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With the help of a multi-configurational Green's function approach we simulate single-electron 
Coulomb charging effects in gated ultimately scaled nanostructures which are beyond the scope 
of a selfconsistent mean-field description. From the simulated Coulomb-blockade characteristics 
we derive effective system capacitances and demonstrate how quantum confinement effects give 
rise to corrections. Such deviations are crucial for the interpretation of experimentally determined 
capacitances and the extraction of application-relevant system parameters. 
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One of the major challenges for the simulation of elec- 
tronic transport in low-dimensional nanostructures con- 
sists in an adequate description of the Coulomb inter- 
action. In ultimately scaled semiconductor nanosystems, 
only a few electrons constitute the current. Thus, the de- 
tails of the electron-electron interaction may become cru- 
cial for the electronic properties, as can be seen in exper- 
imentally observed single-electron chargingeffects in car- 
bon nanotubes and III/ V nanowhiskers [l|, |2| • The chan- 
nel region of a typical application-relevant nanostructure 
involves on the order of 100 single-particle states (or 
sites) and can be highly inhomogeneous and anisotropic. 
Here, external contacts and gate electrodes in general in- 
troduce non-linear perturbations to the system. Hence, 
highly idealized interacting few-level models with a small 
number of effective parameters (constants) become inad- 
equate. For a realistic simulation of quantum transport 
in application-relevant nanodevices, various approaches 
have been used, differing in the formulation of Coulomb 
interaction and contact coupling. While the orthodox 
theory |3, Ij] correctly describes few-electron effects such 
as single-electron tunneling in quasi-isolated quantum 
dot systems, the nonequilibrium Green's function formal- 
ism (NEGF) [a, la, [TJ, Ig also accounts for contact renor- 
malization and dissipation terms. However, the NEGF 
becomes unable to describe the details of charging effects 
of a few fluctuating electrons as soon as a mean-field ap- 
proximation is introduced. 

Recently, we have presented a multi-configurational 
approach (MCSCG) [2|, which represents a consistent 
extension of the mean-field NEGF for the inclusion of 
single-electron effects under application-relevant condi- 
tions. This approach combines the simplicity and scala- 
bility of the mean-field NEGF with a many-body Fock 
space description of the Coulomb interaction of those 
electrons that are resonantly trapped within the nanos- 
tructure. As a key element, the algorithm identifies "rel- 
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FIG. 1: Schematic sketch of the considered coaxially gated 
nanosystem (see text for more details). 



evant" trapped single-particle states that are subject to 
occupation fiuctuations. Within the resulting relevant 
Fock subspace, many-body "configurations" k and asso- 
ciated weights Wk are defined as eigenstates and eigenval- 
ues of the projected many-body statistical operator p'^^ , 
respectively. Consequently, the system Green's functions 
are written as weighted averages: G = ^,^WkG[k], where 
G[k] follows from a Dyson equation with a suitable self- 
energy approximation E[k]. As a self-consistency con- 
dition, p^^^ is obtained by use of the projected many- 
body Hamiltonian H^^^ and the subsidiary condition 
that G^^ should resemble G within the relevant sub- 
space, where G'^^^ is calculated via p*^-^ and H^^^ . 

In this article, we show how the MCSCG formalism 
can be used to calculate electronic transport proper- 
ties of gated one-dimensional (ID) nanosystems in the 
single-electron tunneling regime beyond the mean-field 
approach. Furthermore, we extract realistic system ca- 
pacitances from the simulated TV characteristics. Three 
length scales are relevant for the Coulomb interaction 
energies of the system: The geometrical length L, the 
electrostatic screening length A, and the de Broglie wave- 
length of the electron. The latter becomes important in a 
nanoscale system with quantum confinement, giving rise 
to significant deviations from the expected capacitances. 



In the following, we discuss a ID nanostructure (i.e., 
one lateral mode only) as depicted in Fig.^ In the shown 
coaxial geometry, the channel is surrounded by a metal- 
lic gate electrode, separated by an insulating oxide layer. 
The ID channel of the nanostructure is connected to two 
Schottky contacts (denoted as source and drain) which 
serve as electron reservoirs |l5l[ . In order to visualize the 
electronic structure of the system, a radial cut through 
the inner region of the structure has been inserted, show- 
ing the groundstate probability density |'!/'(2;)|^ and the 
associated local density of states (LDOS) for 1 electron 
inside the channel. For the discussion of interaction ef- 
fects we will now consider an example of an ultimately 
scaled semiconducting nanocolumn (see also [2|]) with an 
effective channel length oi L ~ 19nm, a channel diame- 
ter of dch — 4nm and an oxide thickness of dox = lOnm. 
Furthermore, we chose ech = 15 and eox = 3.9, yielding 
a screening length of A « 3.7nm. The latter is obtained 
from 
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for a coaxial system layout [lOj . The Schottky source and 
drain contacts exhibit a barrier of ^b = 0.5eV (with re- 
spect to the Fermi energy) , whereas the effective electron 
mass is taken as m* = O.OSme. This parameter set was 
chosen in order to clearly visualize the Coulomb effects, 
at the same time being close enough to experimental re- 
alizations. 

Fig.|3a) shows a simulated grayscale map of the drain 
current as a function of the applied gate and drain volt- 
ages in the single-electron tunneling regime (T = 77K) 
[ig. One can clearly identify the typical diamond- 
like shaped black areas of Coulomb blockade with a 
suppressed drain current. Dark-gray areas correspond 
to electronic transport through the lowest longitudinal 
resonantly-trapped state, whereas brighter areas involve 
the first (and higher) excited resonances. Such a be- 
havior is beyond the capabilities of a mean-field NEGF 
that is based on a selfconsistent potential, demonstrat- 
ing the strengths of the employed multi-configurational 
approach. In order to visualize the interaction mecha- 
nism that goes beyond mean-field. Fig. [2Ib) shows the 
local density of states (LDOS) within the channel region 
at different voltage conditions and the associated mean- 
field potentials. For the empty channel "0", we obtain a 
set of quasi-bound states and resulting resonances. With 
one electron inside the channel near equilibrium point 
"1", one electron occupies the lowest state and all fur- 
ther states are shifted to higher energies by the Coulomb 
interaction. Note that the electron does not experience 
its own Hartree potential. Finally, under strong non- 
equilibrium bias conditions at point " X" , the spectral 
properties result from a weighted mixture of an empty 
channel and a channel containing one electron. Here, 
the electronic properties clearly deviate from a mean-field 
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FIG. 2: (a) Simulated single-electron tunneling characteristics as 
a function of the applied gate and drain voltages for T =77K. Here, 
the absolute value of the drain current is mapped to a non-linear 
grayscale. Within the black diamonds of Coulomb blockade the 
electron number becomes quantized as indicated by white numbers. 
AVfj and AVd denote the width and half height of the first dia- 
mond, respectively, (b) Corresponding local density of states and 
mean-field potential within the channel for selected points 0,1, X 
in map (a). Point X corresponds to a nonequilibrium state with 
electron number f» 0.3. 



picture. 

For the considered example, the height of the Coulomb 
diamond with 1 electron is determined by the Coulomb 
interaction of two electrons with opposite spin within the 
lowest channel resonance (groundstate). In contrast, for 
the second diamond one has to consider not only the in- 
teraction between the lowest and the first excited states 
but also the difference between the single-particle ener- 
gies of the two involved states. Hence, its size is signifi- 
cantly enlarged. These observations strongly depend on 
the degeneracy of the considered channel resonances and 
on the ratio of quantization and Coulomb energies. In 
ultimately scaled systems that are addressed in this let- 
ter, the longitudinal quantization energy dominates and 
it is justified to consider the groundstate Coulomb ma- 
trix element only. In contrast, for longer channels, the 
quantization energy may become much smaller than the 
charging energy and we could end up with a situation 



where the occupation of higher excited states (with a re- 
duced Coulomb matrix element) and a re-organization 
of electrons might become energetically favorable for the 
addition of the next electron pjl (and Wigner molecules 
might form [l^)- In the following, we restrict ourselves 
to the discussion of the lowest Coulomb diamond which 
corresponds to the onset of the current flow. 

In general, effective quantities such as the gate capac- 
itance Cq and the total system capacitance C^ can be 
extracted from the shape of the observed Coulomb di- 
amonds [13| (provided that the latter is solely due to 
Coulomb interaction): From the total width AVg and 
half height AVd one obtains Cq = e/AVc and C^ ~ 
e/AVo, respectively. In our example (Fig. ISJ, we have 
AVg « 114mV and 2 AVd « 186mV, and therefore C£ w 
1.41aF and C^ w 1.72aF. Consequently, a Coulomb in- 
teraction energy of E'^ = e^/C^ « 93meV can be derived 
(which is identical to a two-particle Coulomb matrix- 
element, see Ref. |2|). Note that for such ultra-short 
nanosystems, where the quantization effects play an im- 
portant role, these capacitances in general are voltage 
dependent quantities, depending on the wave functions 
of the involved longitudinal channel resonances that are 
being filled at the considered electrode voltages. 

For the classical capacitance of the considered coaxial 
geometry of the ID channel region one obtains 
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Employing the Coulomb Green's function (for details see 
Ref. (sl) for fixed potential boundary conditions, the 
Coulomb interaction energy S^°™ of a homogeneously 
distributed test charge e with an identical background 



charge thus reads as E^° 
homogenous system capacitance 



7C^°" with the total 
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(Note that C^"™ > Cox and C^°" -^ Cox for L/\ -^ oo.) 
Using the given parameters, we thus obtain Cox = 2.30aF 
and Cj^"™ = 3.75 aF. In this context, we therefore define 
the quantity 
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which we refer to as the "quantum confinement factor" 
for the given geometrical structure and parameters (volt- 
ages, etc.). In the example above, we have /e ~ 2.2. The 
observation that /s becomes > 1 for ultimately scaled 
structures can easily be understood if one takes a look at 
the quasi-bound single-electron wave function (Fig.l^b)) 
of the considered Coulomb diamond: Typically, this wave 
function exhibits an enhanced probability amplitude in 



the center region of the ID channel due to quantum con- 
finement in transport direction which is a result of the 
Schottky-barriers at the source and drain contacts. In 
other words, the electron is spatially confined within a 
reduced effective channel length < L. This implies an 
increased Coulomb interaction energy E^ = 93meV com- 
pared to E^°"^ — 43meV for a homogenous charge distri- 
bution, or equivalently C^ < C^°"^ [l3|. Consequently, 
larger voltage changes are required to add one further 
electron to the channel region. This implies that care 
has to be taken with the extraction of geometrical pa- 
rameters from experimental data as will be discussed in 
more detail below. 
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FIG. 3: (a) Simulated capacitances for the first two electrons as a 
function of the channel length L with an adaptive gate voltage. The 
screening length A a; 3.7nm and the eff'ective mass m* = O.OSrrie are 
kept constant, (b) Quantum confinement factor /j] = C|^°'"/C|, 
and gate efficiency rio = ^gI^t, ^^ derived from (a). 

Fig. Ota) shows the simulated capacitances for the 
first two electrons (groundstate wave function and next 
higher state) as a function of the channel length L for 
a fixed A and m* . Here, the gate voltage Vq is cho- 
sen adaptively such that the groundstate energy coin- 
cides with the Fermi energy of the source contact (for 
a constant Schottky barrier height $b = 0.5eV with 
Vd = Vs = 0). For each L, the wave function and the 



associated Coulomb interaction energy E^ = e^ /C^ have 
been obtained via numerical diagonalization of the chan- 
nel Hamiltonian [9j. As expected, for larger L/\ we ob- 



serve that C? 



C*Q and 0^°"^ 



Cox ) which means that 



the source and drain influence becomes negligible in this 
limit. On the other hand, for the ultra-short structure 
with L < 20nm, one can clearly identify the increasing 
influence of the source and drain capacitances. There- 
fore, an increased gate voltage becomes necessary to pull 
the lowest quasi-bound state down to the Fermi energy 
of the source contact for a fixed Schottky barrier height 
(e.g., Vg = 0.52V for L = 50nm, and Vg = I.IV for 
L = lOnm). This is described quantitatively by the effi- 
ciency rjG = Cq/C^ = AVd/AVg of the gate electrode 
at changing the considered quasi-bound energy level, as 
plotted in Fig. Elb). 

Most noticeably for all ranges of L/X, there is a sig- 
nificant deviation of C^ from C|5°™ (and Cq from Cox) 
as can be seen in Fig. El^a). This observation, in fact, 
is described by the quantum confinement factor /s as 
plotted in Fig. |3^b). In the short-channel regime, we ob- 
tain /s > 1, which can be explained by the longitudinal 
confinement of the groundstate wave function and the 
resulting increased Coulomb interaction. Furthermore, 
as can be seen in Fig. |31 the onset of the short-channel 
regime, where rjG deviates significantly from 1 due to in- 
creasing source and drain influence, is shifted to smaller 
channel lengths L < 30nm compared to L > 50nm where 
(jhom gig^j-^g ^Q deviate from Cox- Again, this finding 
can be attributed to the longitudinal confinement of the 
wave function with a reduced electron charge density at 
the source and drain contacts. For larger channel lengths 
L, however, we observe a rapid decrease of /s (or equiv- 
alently, a strong increase of C^ and Cq). This finding 
can be attributed to the onset of a Wigner regime due to 
a significantly reduced average electron density (because 
of an increasing spatial separation of the two considered 
electrons compared to A). 

The discussed confinement effect becomes very impor- 
tant for the experimental determination of application- 
relevant capacitances from measured Coulomb diamonds. 
Analogous to the interpretation of our simulation results 
above, effective capacitances C^ and Cq can be calcu- 
lated from the experimentally determined values of AVd 
and AVg at low temperatures (typically AksT < E^). 
One has to be careful with the interpretation of these 
effective capacitances as far as the derivation of geomet- 
rical parameters from the latter and their relevance at 
room-temperature are concerned. For example, it is pos- 
sible to calculate dox from C^. By inspecting Eqs. 1^ 
0}, however, one can see that this is a non-trivial task 
since A depends on dox, and /s (which results from a 
Coulomb matrix-element) in turn depends on A. If we 
identified C^°" with C^ = e/AVo (i.e., setting /s = 1) 



from experimental data, we would overestimate dox- In 
the considered example, we would end up with dox > 
50nm (instead of lOnm) which is a consequence of the 
logarithm in Eqs. I|1I2|I . 

In conclusion, we have demonstrated that the multi- 
configurational selfconsistent Green's function approach 
(MCSCG) is able to model single-electron Coulomb 
charging effects in ultimetely scaled gated ID nanostruc- 
tures, providing a means for the determination of inter- 
action energies and resulting effective capacitances. In 
turn, a quantitative discussion of quantum confinement 
corrections to the system capacitance was given, resulting 
in a significant deviation from the classical capacitance 
which has an immediate consequence for the interpreta- 
tion of experimental system parameters. 
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In this example, up to 64 configurations of relevant states 
with a diagonal interaction approximation were consid- 
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The reduction of C^ with respect to C^"™ is not to be 
confused with the concept of a finite quantum capaci- 
tance Cq = e^D jEp ) with a thermally averaged density 
of states D(Ef) uM- The considered nanosystem in the 
single-electron regime exhibits a quasi-zero-dimensional 
electronic structure, where Cq equals or oo. 



